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ANALYSIS OF LONGITUDINAL PILOT-INDUCED 


OSCILLATION TENDENCIES OF YF-12 AIRCRAFT 

John W . Smith and Donald T . Berry 
Flight Research Center 


INTRODUCTION 


Ever since the first manually controlled flight of powered aircraft , pilots have 
experienced pilot- induced oscillations (PIO) (refs. 1 and 2) . Although their cause 
and nature differ from one aircraft to another, the problem is usually caused 
by the pilot's phasing his control inputs so that the aircraft's aerodynamic or con- 
trol modes, or both, are amplified instead of attenuated. Two examples of longitu- 
dinal PIO's encountered with earlier production aircraft are force feedback directly 
to the pilot due to a bobweight (ref. 3) and force feedback through the cables 
and linkages on a reversible system due to control surface balances (ref. 4) . More 
recently, nonlinear elements in stability augmentation systems (SAS) have aggra- 
vated the pilot/control system interaction problem by adding more lag, which can 
contribute to the pilot's initiating and sustaining an oscillation (ref. 5) . 

Despite the long history of PIO's , a generally valid prediction technique does 
not exist. One of the means available to designers for their prediction and preven- 
tion is to compare the characteristics of new aircraft with the characteristics of air- 
craft that have proved to be resistant or prone to these oscillations . 

To contribute to the body of knowledge on aircraft PIO characteristics , this 
report analyzes PIO's experienced with YF-12 aircraft. Two types of PIO's have 
been experienced. The first, which is more common, takes place during aerial re- 
fueling when the pilot gain is high. The pilot tends to interact with the aircraft's 
SAS-on short-period and structural modes . This naturally makes the task more 
difficult, even though the attitude disturbances are not large. The second type of 
PIO discussed, although more infrequent, is considerably more severe in terms of 
aircraft rates and accelerations. Again, the PIO occurred during refueling, where 
pilot gain was high . 

This report discusses aircraft ground and flight tests, analog simulator studies, 
and theoretical analyses that were made at the NASA Flight Research Center to study 
the phenomenon . 



SYMBOLS AND ABBREVIATIONS 


Physical quantities in this report are given in the International System of 
Units (SI) and parenthetically in U.S. Customary Units. The measurements were 
taken in U.S. Customary Units. Factors relating the two systems are presented in 
reference 6. 
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frequency , cps , or function 
describing function 

mechanical gearing between stick and elevator 

gravitational constant , m/sec^ (ft/sec^) 
total hysteresis width 
Laplace transfer function , feedback loop 
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mass , kg (slugs) 
nonlinear element 
dynamic pressure , N/m^ (Ib/ft^) 
reference wing area, m^ (ft^) 

Laplace transform 
transfer function 
time , sec 

velocity, m/sec (ft/ sec) 

output quantity 

angle of attack , deg or rad 

flightpath angle , deg or rad 

rate of change of flightpath angle, rad/sec 

incremental change or characteristic equation 

elevator deflection, deg or rad 

pilot elevator deflection command, deg or rad 
stability augmentation system input , deg or rad 
stick deflection, deg or rad 


^ damping ratio 

0 pitch attitude , deg or rad 
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pitch attitude error , ~ 

pitch rate , deg/sec or rad/sec 
pitching acceleration, deg/sec^ or rad/sec^ 
real part of a root, rad/sec, or root mean square 
time constant 

phase angle , deg or rad 
frequency, rad/sec 

critical 

center of gravity 
closed loop 
crossover 
cockpit location 
initial 
inboard 

natural frequency 
output 
outboard 
reference 

sequential calculation in the order indicated 


VEHICLE DESCRIPTION 


Airplane 

The YF-12 airplane (fig. 1) is an advanced, twin-engine, delta-wing interceptor 
designed for long-range cruise at Mach numbers greater than 3 and altitudes above 
24,400 meters (80,000 feet) . Airplane physical characteristics are given in refer- 
ence 7 . 
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Figure 1. Three-view drawing ofYF-12 airplane. 


Flight Control Systems 

Pilot control inputs are transferred by cables and linkages to servos in the aft 
section of the airplane . The inputs from the pilot and the SAS are summed in the 
servos that control the inboard elevon power actuators. The position of the inboard 
elevons is transferred by rods and linkages across the nacelles and used as an input 
to the servos that control the outboard elevon power actuators . 

Limited-authority rate dampers are utilized in all three axes to provide the basic 
airplane with additional stability and damping. Damper position limits in pitch are 
2.5° trailing edge up and 6.5° trailing edge down. The maximum damper surface 
rate is 15 degrees per second . In the feedback loop , the SAS employs gain sched- 
uling that varies according to altitude and dynamic pressure . However, at the flight 
conditions of interest the gain is fixed in pitch at a value of 0.75 deg/deg/sec . 


6 



INSTRUMENTATION 


Standard NASA instrumentation was used to measure the pertinent flight quan- 
tities and control positions . The sensors used to detect accelerations and angular 
rates were near fuselage station (F.S,) 12.7 meters (500 inches) . The data were 
sampled at 200 times per second, conditioned, and recorded both on board and on 
the ground. The data were transmitted to ground stations through a pulse code 
modulation (PCM) system . The PCM system is considered to be accurate within 
2 percent of a full-scale deflection. 


ANALYSIS AND DISCUSSION 


Pilot-Induced Oscillation Experiences 

Aerial refueling is necessary to accomplish the design missions of the YF-12 
aircraft . During the refueling maneuver the pilot must fly the aircraft to a position 
close to a tanker and maintain that position for several minutes . The task is 
demanding and requires tight , attentive attitude control on the part of the pilot . In 
general, as compared with other aircraft, the handling qualities of the YF-12 air- 
craft during refueling are considered to be good . Nevertheless , PIO tendencies 
during refueling have been reported. As described by the pilot, the reason for 
this tendency is a hobbling motion at the cockpit , which is caused by the flexing of 
the aircraft . The hobbling motion is benign , in that the oscillations are of small 
amplitude . However , the attitude changes act as triggering cues under conditions 
of high pilot gain and result in a coupled interaction . The interaction tendency 
increases as more and more fuel is taken on board the aircraft. 


Figure 2 is a time history of the 
latter portion of a refueling maneuver . 
During this maneuver the control input , 
6 , was less than 3° peak to peak at a 

s 

frequency of approximately 1 cycle per 
second. These control demands, which 
are somewhat annoying, are typical of 
the control inputs that the pilot has to 
make to keep the pitch excursions to a 
minimum . Although small , the control 
inputs are continuous , particularly dur- 
ing the final portion of the refueling 
maneuver , and they require additional 
effort on the part of the pilot. 

The second type of PIO experienced 
is associated with large control inputs 
and large -amplitude responses. This 
type of PIO has occurred at least twice . 
In one instance, the PIO was triggered 
by an overshoot in longitudinal trim 
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Figure 2. Time history during refueling. 
Typical of small-amplitude PIO tendency. 
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because of a faulty trim switch. This happened as the aircraft was approaching the 
tanker just before hookup. The pilot reacted normally to the trim overrun and took 
abrupt corrective action to keep the airplane from reaching its g limit . This caused 
a pitch oscillation that persisted for several cycles (fig. 3) . In this particular case, 
the airplane's instrumentation was limited, and the only variables available are 
shown. It is apparent, however, that excursions of approximately -Ig to 3g 
occurred in normal acceleration . Control surface deflections were not recorded , 
but the pilot reported the incident as a definite PIO . 

A study was initiated to determine the basic cause of these control difficulties . 




■2 * ' 1 I I 

0 12 3 4 


t, sec 


Figure 3. Large -amplitude PIO time history. 


Small-Amplitude Pilot-Induced Oscillations 

The PIO's experienced during aerial refueling (fig. 2) were of relatively small 
amplitude. Although small, such oscillations can be bothersome, because once they 
are excited by control activity, additional control inputs may be required on the 
part of the pilot to minimize them. The oscillations resulted from the coupled dynam- 
ics of both the rigid body motion and the flexing of the aircraft , and the pilot found 
it difficult to distinguish between the modes . 
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Because of the small amplitudes involved, linear analysis was expected to be ade- 
quate to analyze the problem. The model shown in figure 4, which shows the elements 
of the longitudinal control system and their relationship to the pilot, was used for the 
analysis . The first structural bending mode was included in the model because of 
the airplane's tendency to bobble and flex. In the figure , 0^ refers to the pitch angle 

sensed at the pilot's station, whereas 0 is the rigid body pitch angle. The SAS reacts 


First bending mode 



Figure 4. Linear model descriptive of longitudinal control system. 

to what is essentially the rigid body pitch angle, since the system's gyroscope is on 
a zero slope point of the bending mode deflection curve (fig. 5, based on ref. 8) . 


Fuselage station, In. 

0 2 4 6 8 10 12 14x10^ 


Normalized 

deflection 



Figure 5. First structural mode shape. 
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Consequently, the contribution of the bending mode to pitch angle can be repre- 
sented as summing downstream of the SAS feedback loop , as in figure 4 . The con- 
trol system's characteristics and the aerodynamics used in the analysis are de- 
scribed in appendixes A to D. The poles and zeros that represent the closed-loop 
or equivalent pilot-vehicle system are analyzed in linear fashion by the methods 
used in reference 1 and are illustrated by the root locus in figure 6. Increasing 
gain, K„, which is the same as increasing pilot gain, drives the poles of both the 



a, rad/sec 

Figure 6. Root locus of longitudinal system indicating effect of increasing 


0 

gain in the attitude loop. ^ _ 

e l+A~sH(s) 

P O 

e 

first bending and aircraft SAS-on short-period dynamics to a more oscillatory 
condition. In fact, if the gain is large enough, the short-period mode would tend 
toward instability near 7.8 radians per second. This depicts the less severe 
small- amplitude PIO's shown in figure 2. The interaction of the pilot with the com- 
bination of the airplane's SAS-on short-period mode and the structural first bend- 
ing mode zeros is responsible for the added control effort during the refueling 
maneuver. The root locus also indicates that the first bending mode goes unstable 
at a lower pilot gain than the short-period mode . However , the frequency of the 
first bending mode is far enough removed from the frequency of the pilot's control 
input so that this does not happen. This was verified by random analyses like 
those in figures 7 and 8 of data recorded while the airplane was connected to and 
disconnected from the tanker . 
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Figures 7 and 8 show values of amplitude spectral density per cycle-per-second 
filter bandwidth over the pertinent frequency range . This method of analysis is 
discussed in reference 9. Figure 7 compares stick deflection versus frequency 
during refueling with stick deflection when the airplane is controlled just enough 
to maintain altitude. A comparison of the two curves shows that the pilot gain is 


Connected 



Figure 7. Comparison of amplitude spectral density variation of primary pitch 
control input during flight connected to and disconnected from the tanker. 


much higher during refueling than when he is just maintaining altitude, even 
though the root-mean- square variation in pitch attitude is approximately the same 
(fig. 8) . Furthermore, despite the peak in pilot input at 1 cycle per second, which 
is attributable to the PIO tendency , the pilot input falls off sharply at frequencies 
higher than 1 cycle per second (fig. 7) . The 1-cycle-per-second frequency is less 
than one-half the frequency for the structural first bending mode (fig. 6); con- 
sequently , the pilot is unlikely to excite or sustain pure bending mode oscillations . 
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Og, deg/cps 



Figure 8. Comparison of amplitude spectral density variation of pitch attitude 
during flight connected to and disconnected from the tanker. 


Large-Amplitude Pilot-Induced Oscillations 

As shown in figure 3 , during the large-amplitude PIO the aircraft experienced 
excursions of -Ig to 3g at a frequency of just over 0.5 cycle per second. This 
frequency is lower than predicted by linear analysis , and the magnitude of the 
responses suggests that the system's nonlinearities must be included in an analysis 
of this PIO tendency . 

Although there are no control input data for the incident shown in figure 3, sim- 
ulator tests , airplane closed-loop ground tests , and limited flight tests confirmed 
that the SAS limits were exceeded. The tests also verified the system's characteris- 
tics during PIO's (appendix A) . Figure 9 shows a model of the airplane plus SAS 
that was derived from information obtained from the tests . The nonlinearity of the 
description was necessary for the analysis because of the rate and position limiting 
that occur in the SAS loop during a severe PIO. A digital computer program was 

0 \ 

written to give the closed-loop frequency responses of g — and g — at both the cen- 

e e 

P P 

ter of gravity and the cockpit , in addition to minimum boundary curves , regions , 
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and intersects for nonlinear operation. Linear conditions, which were used to make 
checks and comparisons , were also programed on the digital computer . (An outline 
of the digital computer program is presented in appendix E, along with the transfer 
functions used to describe the system's elements.) 



Figure 9. Model used to describe nonlinear airplane 
pilot open-loop longitudinal response. 


Figures 10(a) and 10(b) show the results of pitch attitude response per com- 
mand input versus frequency. Three calculations are presented. One is for a 
linear system, and the other two are for a nonlinear system. The calculations for 
the nonlinear system are for a control input ,8 , of an amplitude typical of a PIO 

and for different damper position limits. The SAS rate limit is 12.6 degrees per 
second. Amplitude attenuation occurs in a normal fashion (fig. 10(a)) . A compari- 
son of the systems shows that the normal attenuation is reduced by nonlinearities , 
particularly near the PIO frequency, up to approximately 4 radians per second. A 
sizable increase in the position limits (from 2.5° to 6.5°) makes little difference in 
the closed-loop gain at the PIO frequency , indicating that rate limiting is the dominant 
influence. Above 5 radians per second, nonlinearities have virtually no effect. (The 
peak at 15.7 radians per second is due to structural first bending . ) 


0 

Figure 10(b) shows the corresponding phase angle variation of g — versus 

e 

P 

frequency for the same three conditions. Again, nonlinearities add to the phase 
lag in the vicinity of the PIO frequency , and above 5 radians per second there is 
less difference between the three systems . The increase in phase lag and lower 
attenuation in the PIO frequency range would increase the tendency for PIO's. 


13 








(a) Amplitude versus frequency . 


Linear 

Position limit = 2.5° I 

Position limit = 6.5° ) 


Nonlinear 



(b) Phase angle versus frequency . 


Figure 10. Comparison of linear and nonlinear calculation of aircraft closed- 
loop pitch attitude near PIO conditions . Rate limit = 12. 6 degrees per second 
at the cockpit; pilot loop open; 8 ^0.1 radian (half the peak-to-peak ampli- 
tude) for nonlinear cases . p 
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Presented in figure 11 is the incremental change in normal acceleration at the cen- 
ter of gravity per commanded input versus frequency . The three calculations pre- 
sented are those shown in figure 10 . One is for a linear system and two are for a non- 
linear system. The calculations for the nonlinear system are for different damper 
position limits and a control input, 6 , of an amplitude typical of a PIO . The SAS 

rate limit is 12.6 degrees per second. The results of the calculations are compared 


o Airplane ground tests 
D Simulator studies 
O Flight tests 



0, rad/sec 

Figure 11. Comparison of experimental data with calculations of normal-force- 
sensitive parameter versus frequency. At the center of gravity with no body 
bending; SAS rate limit - 12. 6 degrees per second; 5^ ^0.1 radian (half the 

peak-to-peak amplitude) for nonlinear cases. p 


with the results of airplane ground tests, simulator studies , and flight tests . The in- 
flight PIO's occurred at the conditions listed in table 1. The test data for the larger 
input amplitudes (solid symbols) , most of which were obtained from ground tests 
and simulator studies for safety reasons, tend to peak out near a frequency of 
3 radians per second. The calculated data predict a somewhat lower PIO frequency 
than the experimental data, probably because the experimental data were obtained 
for slightly different flight conditions than were assumed for the calculations . Sig- 

Aa 

nificantly , both the calculations and the experimental data show a rise in — near 

s 

the frequencies where PIO's were experienced, which means that the airplane is 
more responsive to control inputs in this frequency range . Aside from showing the 
tendency for PIO's, the comparison of the experimental data with the calculations 
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shows that for large input amplitudes the closed-loop nonlinear model describes 
the actual input-output relationship well. It also shows that at low frequencies the 
magnitude of position limiting has a sizable effect on the input-output relationship . 
At the PIO frequency (3.14 radians per second) , however, the difference between 
the linear and nonlinear models is due primarily to rate limiting . 

TABLE 1.— MAXIMUM PEAK-TO-PEAK AIRCRAFT RESPONSE AND CONTROL 
INPUTS DURING PILOT-INDUCED OSCILLATIONS 

(a) Flight tests 



Pilot-induced oscillation 

First half cycle 

Completely developed 

A cps 

0.75 

0.61 

Aa , g 

2 (0.2 to 2.1) 

4 (-1.0 to 3.0) 

A0 , deg 

4.9 

10.4 

A9 , deg/sec 

13.0 

26.5 

A6, deg/sec^ 

71 

95 


(b) Ground tests and simulator studies 



Pilot-induced oscillation 

First half cycle 

Completely developed 

Ground 

tests 

Simulator 

studies 

Ground 

tests 

Simulator 

studies 

A6^ , deg 

11.5 

10.6 

13.4 

10.4 

A6 , deg 

13.8 

16 

18.6 

16 

f. cps 

0.80 

0.80 

0.61 

0.56 

(v) 

^ P^SAS off 

1.47 

1.46 

1.66 

1.46 

A6^ , deg 
P 

9.5 

11 

11.2 

11 
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The airplane closed-loop boundaries between linear and nonlinear operation, 

S = /■(N- , , co) , were determined by using the transfer functions descriptive of the 

hardware and the nonlinearities for the SAS . Figures 12(a) and 12(b) present the 
minimum conditions for position limits of 2.5° and 6.5°, respectively, for rate limit- 
ing of 15 degrees per second . The level of the average control input in these figures 
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(a) Position limit = 2 . 5° . 



(b) Position limit = 6.5°. 


Figure 12. Limiting envelopes as a function of input amplitude and frequency . 
Rate limit - 15 degrees per second. 
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was obtained from records that were taken continuously during refueling and ana- 
lyzed spectrally, as previously discussed. The analysis showed that most of the 
time the pilot's control inputs were below the limiting boundaries. In fact, if it is 
assumed that there is a Rayleigh distri- 
bution for 6 , even at the highest aver- 


age input the probability that a control 
input will not exceed 0 . 05 radian is on 
the order of 95 percent. Therefore, it 
is highly unlikely that the pilot will get 
into a rate- or position-limiting situa- 
tion. However, during the latter por- 
tion of the refueling operation (at a 
heavy condition) the crew has com- 
plained about a tendency toward PIO's, 
and on one occasion , although a PIO 
did not develop fvilly, the pilot's con- 
trol input, 6 , exceeded 7° (t = 0 to 

s 

2 sec, fig. 13) or 


8 =0.07 radian at a 


frequency of approximately 3 radians 
per second for two or more cycles . 
Figure 12 also shows the conditions 
under which a PIO did develop fully . 
The peak-to-peak values of control in- 
puts and aircraft responses are listed 
in table 1 . In general , for the fre- 
quency range of interest, rate limiting 
is more prevalent and more of a prob- 
lem than position limiting. 

The calculations from the digital 
computer program in figures 10, 11, 
and 12 are typical and are felt to be 
descriptive of the characteristics of the 
airframe and the augmentation system . 
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Figure 13. Time history during re- 
fueling. Typical of early part of 
severe PIO. 


Closed-Loop Nonlinear Analysis With Pilot in the Loop 

Figure 14 presents a nonlinear model of the airplane with the pilot in the loop . 
The nonlinearities of the artificial feel system are included in the model , as well 
as the nonlinearities of the SAS . As described in reference 1 , a closed-loop analy- 
sis can be performed by considering the pilot to be controlling to some desired 
attitude, which minimizes pitch attitude error, 0^. To sustain an oscillation in a 

closed-loop system, the total loop gain must be equal to 1 at a phase shift of -180°. 
If the pilot's adaptation during a PIO is assumed to be synchronous (in phase) with 

F 

g 

the pitch attitude error signal, ^ can be considered a pure gain (ref. 3) . 

e 
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Figure 14. Nonlinear model of pilot-aircraft closed-loop system. 


Therefore , for the PIO analysis , 

^ (Describing function) ^(A/C)^^J=-1 

where the describing function is a combination of the force hysteresis loop plus the 
nonlinear gearing characteristics (appendix F) . Both the describing function and the 
transfer function of the airplane plus the SAS , (A/C)^p are a function of both the 

frequency and initial amplitude of the pilot input , 6 ; therefore , the pilot gain re- 

quired for PIO , -g — , is also a function of 6 . Figure 15 shows the variation of PIO 

e 

frequency and pilot gain as a function of input amplitude for a -180° phase shift and 
a total gain of 1. The PIO frequency is reduced to 50 percent to 75 percent of the 
linear crossover frequency (approximately 7.8 radians per second) . 

The conditions that triggered PIO's in flight are shown in figure 15 and described 
in table 1 . The trend of the flight data in figure 15 is similar to the trend of the non- 
linear calculations , but the PIO occurred at a slightly lower frequency than pre- 
dicted . As shown in the lower plot , at the PIO input amplitude (6 =0.1 rad) , the 

pilot gain necessary to sustain a PIO is only half the gain required to sustain a PIO 
with a linear system . Therefore , if the pilot gain was high , as during formation 
flight, which minimizes 0^, a large enough control input could indeed trigger a PIO. 

Such an incident as runaway or overrun trim could easily result in such an input. 

It was felt that increasing the SAS rate limit to that of the servoactuator (from 
12.6 degrees per second to 30 degrees per second) might alleviate the problem. 
Consequently , calculations were made to see what could be expected with the higher 
rate limit. For the condition calculated, which was near the PIO input amplitude, 
the pilot gain necessary to trigger a PIO with the nonlinear system was only 30 per- 
cent lower than the gain that would trigger one with a linear system . In addition , 
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with the higher rate limit, the PIO frequency increased to approximately 6.50 radians 
per second. These results imply that a SAS loop that was more responsive at the 
larger amplitudes would reduce the PIO tendency . 



3x10^ 


2 


F^/e 


e' 


Ib/rad 
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0 


Figure 15. PIO frequency and minimum pilot gain required for PIO versus 


input amplitude for 



= 0. 000944 rad/N (0. 0042 rad/lb) . 


CONCLUDING REMARKS 


Two types of pilot-induced oscillation (PIO) conditions were encountered dur- 
ing flight tests on the YF-12 airplane. To analyze the problem, studies were con- 
ducted using an analog simulator , ground tests , and linear and nonlinear calcula- 
tions . 

A small- amplitude PIO tendency near a frequency of 1 cycle per second existed 
when the pilot gain was high during refueling. The demanding nature of the control 
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task and the response characteristics of the airplane (short-period poles) , coupled 
with structural bending' (zeros) , resulted in a continuous but small control input 
type of PIO . The structural flexing of the airplane , predominantly the first bending 
mode , was apparent to the pilot and caused him some concern at first after a dis- 
turbance; however, he avoided controlling continuously in or out of phase with the 
first bending mode response (2 . 5 cycles per second) . During the early portion of 
refueling or at the lightweight condition , this PIO was only a minor annoyance . 

A large- amplitude (-Ig to 3g) and potentially more serious PIO tendency is 
caused by SAS rate and position saturation , which in turn cause phase lag and less 
amplitude attenuation . These nonlinearities can reduce the pilot gain required to 
sustain a PIO by more than 50 percent . SAS rate limits are more detrimental than 
position limits . The pilot normally does not make inputs large enough to cause 
saturation, so this situation is unusual. However, it can be initiated by an un- 
expected disturbance like a faulty trim command. 


Flight Research Center 

National Aeronautics and Space Administration 
Edwards, Calif., December 11, 1974 
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APPENDIX A. -AIRCRAFT MODELING AND VERIFICATION 


GROUND TESTS 


Ground tests were conducted on the YF-12 airplane to determine its eleven, SAS , 
and actuator characteristics . An analog computer was tied in to the SAS to simulate 
the airplane's aerodynamics. Tests were conducted to determine the control inputs 
and SAS nonlinearities during a typical PIO condition, since such testing in flight 
might have endangered the airplane . Briefly , the tests involved applying sinusoidal 
inputs to the augmentation system and control cycles and step pulses to the mechani- 
cal system. 

Figure 16 describes the complete simulation. Airplane systems used in the sim- 
ulation consisted of gearing, G^; actuators and servos, A(s); and feedback gain and 

shaping, H(s) . The inputs to the computer were electrical signals from the trans- 
ducers on the surface actuators . The gains and are the scale factors necessary 

to make the signals proportional to the desired control surface input and gyro output , 
respectively. The computer was programed using a two-degree-of-freedom simula- 
tion. The equations of motion used are presented in appendix D. Consideration of 
the first bending mode was confined to its effect on the normal acceleration at the 
cockpit , which was displayed on a meter in the cockpit . The first bending mode 
dynamics data were obtained from flight tests . The stability coefficients were ob- 
tained from predictions that were adjusted to flight data for a Mach number of 0.77 
and an altitude of 7620 meters (25,000 feet) . Appendix D lists the actual values of 
these coefficients . 


Display 



Figure 16. Description of ground test simulation. 
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APPENDIX A - Continued 


Figure 17 presents typical data from these ground tests. The first part of the 
time history (t = 3 to 7 sec) is a control release. It was used to check the primary 
control response and the SAS gain. From 11 seconds on, the primary control sys- 
tem was pumped at a frequency near the PIO condition (co 4.7 radians per second 
and A6 = 17° peak to peak) . Peak-to-peak excursions of 3.6g occurred in normal 

acceleration. The damper 6 bottomed out at -2.5°, and rate limiting of an 

®SAS 

average value of 12.6 degrees per second was evident throughout most of the cycle. 
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Figure 17. Time history of control and system response 
obtained from ground tests simulating a PIO condition. 
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APPENDIX A - Continued 


CONTROL SYSTEM RESPONSE 


The transfer functions and frequency responses of the inboard and outboard 
elevons to the output of the SAS servo's linear variable differential transformer 
(LVDT) is shown in figures 18(a) and 18(b) . The amplitude ratios are fairly flat out 
to a frequency of 20 radians per second. At that frequency the phase difference 
of the inboard elevon is approximately -25° and that of the outboard elevon is approx- 
imately -100°. The transfer functions, which were calculated by the airplane manu- 
facturer, show good agreement with the data acquired with the actual hardware. 


0,0 Test data 



1 10 100 

u), rad/sec 


(a) Inboard elevons . 


0,0 Test data 



1 10 100 

w, rad/sec 


(b) Outboard elevons . 

Figure 18. Frequency response of the power actuators , 
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APPENDIX A - Continued 


To obtain better short-period closed-loop response, a lag-lead network was 
installed in the SAS feedback loop. The frequency response of this filter, and that 

of its model , - v .— , is shown in figure 19 . The low order gain is approxi- 

mately -2.5 decibels. With increased frequency the model data closely fit the test 
data. Above 12 radians per second, the dynamic characteristics of the servo were 
subtracted from the test data, as shown in the figure. 


o Filter and servo ) ^ ^ j 

a Filter alone ( 

Filter alone, model data, 


0.375(s + 8) 
s +4 



01 , rad/sec 



Figure 19. Frequency response of the lag-lead filter. 


SIMULATOR STUDY 


To acquire further evidence of the PIO problems , an analog computer was pro- 
gramed to simulate the control system and the airplane's aerodynamics. The mech- 
anization of the control system included the SAS nonlinearities position and rate 
limits. The airplane aerodynamics consisted of a three-degree-of-freedom mechan- 
ization. The stability derivatives for the simulation were based on flight test data. 

A time history typical of the results is shown in figure 20. The input to the stick, 

6 , was sinusoidal, with amplitude and frequency as the principal variables. 

In general, the simulation showed that a PIO could be sustained and that the severity 
of the PIO depended on input amplitude and frequency . 

Table 1 summarizes the peak-to-peak variation of the listed parameters obtained 
from ground tests of the airplane and from simulator studies . 


FLIGHT TESTS 


Flight test data of a limited nature were also obtained. The tests consisted of 
the manual cycling of the mechanical controls at various frequencies and at approxi- 
mately one-third the g's experienced during a fully developed PIO. 
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APPENDIX A - Continued 





Figure 20. Analog time history with sinusoidal control input typical of a PIO . 


The time history (fig. 21) indicates that nonlinearities are reached in the 
control system even under these restricted conditions . The data show that the 
damper bottoms out near -2.5° at low frequencies (f ^ 0.50 cycle per second) and 
rate limits at approximately 15 degrees per second at slightly higher frequencies 
(f ^ 1.25 cycles per second) . 
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APPENDIX A ” Concluded 
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0 
-1 


n ^ ^ 
cp 







(a) Damper position limiting, (b) Damper rate limiting, 

f ^ 0,50 cps . f ^ 1.25 cps . 


Figure 21, In-flight time history. Intentional control cycling at refueling flight 
conditions . 
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APPENDIX B. -ACTUATOR CHARACTERISTICS 


The longitudinal control of the YF-12 airplane is provided by four segmented 
eleven surfaces that are faired into the trailing edge of the wing . An output from 
the inboard set of eleven actuators mechanically positions the outboard set. The 
diagram below is representative of the mechanical flow and transfer functions for 
the actuators on either side of the airplane . 




1568 

^ 

+ 50.5s + 1568 








1149 

1 




(s + 33.9)2 



^inbd 


'outbd 


Control power is the product of control effectiveness and surface deflection, as 
follows: 


M = 6 

6 e 
e 

Elevon effectiveness is less on the inboard side than on the outboard side at 
the particular flight condition of interest (M =0.77, h = 7620 m (25,000 ft)): 

C C n p 

m mg 

— = 0.45 and — ^ = o.55. Therefore, Mg is simply 

m„ e 


Mo — Mo + Mo 

o o o 

® ^inbd ^outbd 

The equivalent control power would be given by 

M = M /o.456 + 0.556 \ 

°e\ ^inbd ^outbd) 

where the difference in effectiveness is apportioned to the surface deflection . Sub- 
stituting the transfer functions for 6 and 5 gives 

^inbd outbd 


/ 5 \ 

(0.45) (1568) ^ 

(0.55)(1568) (1149) 

16 ) '^6 
V ®R/ ® 

"r 

s2 + 50.5s + 1568 

(s2 + 50.5s + 1568) (s + 33.9)2 ^ 
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APPENDIX B - Concluded 


By combining both transfer functions , an equivalent fourth-order description re- 
sults: 


^ 705.6 (s^ + 67.8s + 2553.54) 

(s^ + 50.5s + 1568)(s + 33.9)^ 

rC 

In factored form , 

^ 705.6 (s + 33.9 ± j37.474)^ 

(s + 25.25 ±j30. 503)2 (s + 33.9)^ 


(Bl) 


(B2) 
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APPENDIX C. -STRUCTURAL MODE ANALYSIS 


Long, slender aircraft like the YF-12 airplanes exhibit a considerable amount of 
structural bending throughout the flight envelope. Typically, the response of this 
type of airplane to a control input is somewhat modified by the bending that takes place 
while the structural modes are excited. Because of low damping, the longitudinal 
first bending mode superimposes an almost continuous small-amplitude sinusoidal 
response on the overall airplane's responses. For example, after a sudden input to 
the elevens in an upward direction (fig. 22) , the aerodynamic loading on the elevens 
causes the wings to deflect downward. In turn, the mass near the center of gravity 
accelerates downward relative to the rest of the mass along the fuselage, and more 
so at first because of bending caused by the control input . The accelerations due to 
any control action are sufficient to excite the first bending mode . The initial ampli- 
tude of the first bending mode depends primarily on the abruptness and the magni- 
tude of the control input . 



Figure 22. Bending response directions after control input. 
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APPENDIX C - Continued 


The total pitching acceleration , 0 , at any point on the fuselage after a step input 
results from both the airplane's response and the additional bending of the fuselage: 


9 . 


s s +• 


M„ 6 
6 e 
e 


s^6 

B B e 


(Cl) 


+ 2C CO s + (jd^ s^ + 2K„co„s + col 
n n ^B B B 


where K_ is a constant proportional to the additional amount of bending at that loca- 
tion . 

If the initial value theorem is applied to equation (Cl) , the following equation 
results: 


e = M + K cc^ (C2) 

e 

Then the initial value of 0* is measured from a flight record or time history of 9, as 
follows: 



The two slopes are obtained along with the bending mode natural frequency, 
cOg , and bending mode damping ratio, . The value of Mg is usually known from 

e 

wind tunnel or flight test results; therefore, enough information is available to 
estimate K„. 

For example , it was determined from flight data that the first bending mode 
frequency was 15.7 radians per second. It was also determined that 0^ was 1.8 for 

9 9 

the airplane plus the first bending mode and that 0^ was 1. 1 for the airplane alone. 
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APPENDIX C - Concluded 


From earlier flight test data it was determined that Mg was approximately 6.0. 

e 

Making the substitution into equation (Cl) and taking the ratio of the two different 
initial responses, 

^ g 6.00 + Kg246 

TT" Oo 


or 


K 


B. 


= 0.0155 


(Cl) 


instrument location 


For this particular example , the instrument was elsewhere than in the cockpit; 
therefore, the mode shape was needed to determine the value of K at the cockpit. 

Figure 5 presents the normalized bending deflection of the airplane . YF-12 
bending mode data are published in reference 8. All deflections are relative to the 
maximum deflection , which occurs at the nose . The ratio of the slope at the cockpit 
to the slope at the instrumentation package is approximately 1.36, so 


K 


B 


= 1.36K 


cp 


B. 


instrument location 


(C4) 


Therefore , 


K 


B 


cp 


0.021 
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APPENDIX D . -EQUATIONS OF MOTION 


A two-degree-of-freedom representation was considered to be adequate: 


y = L a-^ + L„ 8 
' a V 8 e 

(Dl) 

e 


0 = Mx0 + M a + Mjj 8 
9 a 8 e 

e 

(D2) 

P 

II 

CD 

I 

(D3) 


Taking the Laplace transform and arranging the coefficients in a matrix , 
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from which the transfer functions can be derived as follows; 


9 

6 


L - M Lj, + sMj, 
a 8 a 6 6 

e e e 


e s 


X = 

6 

e 


[s* - (M-L„ M J] 


sA 


The first bending mode is represented by 


(D4) 


(D5) 



^ ^ “b 


(D6) 
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APPENDIX D - Concluded 


The flight conditions and flight test data used for the analysis are listed below: 


g = 9 . 754 m/sec^ (32 ft/sec^) 


hp = 7620 m (25,000 ft) 


= -3.486 per sec^ 


Mg = -6.084 per sec^ 


= 0 .904 per sec 

= -0.609 per sec 

Lg = 0 . 180 per sec 

e V = 239 m/sec (786 ft/sec) 

M =0.77 
n 

From bending mode data and calculations (appendix C) , 


Kg = 0.021 

cOg = 15.7 rad/sec 
^g = 0.05 


For these flight conditions , the transfer functions become approximately 


_e_ ^ -6.00(5 + 0.8) 

s(s^ + 1.50s + 4.00) 

Y ^ 0.180(s^ + 0.60s - 27) 
s(s^ + 1.50s + 4.00) 

and 

0 -5.15 

o 

e + 1.57s + 246 

cp 


(D7) 


(D8) 


(D9) 
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APPENDIX E. -BASIC ANALYSIS AND COMPUTATIONAL PROCEDURE 


The following flow diagram represents the YF-12 airplane's basic control sys- 
tem , aerodynamics , and structural dynamics . 



In this diagram, 6 is the input commanded by the pilot and sums with the 

damper feedback signal, In the SAS loop a lag-lead filter, T^, is implemented 
in addition to the damper rate limits, N^, and position limits, N^- The output param- 
eter, Xg, depends on the transfer function: T^, T^, Tg, and so forth. 

This representation can be used for both linear and nonlinear calculations. For 
linear calculations, and equal 1. For nonlinear calculations, the existing 

amount of rate limiting and position limiting is assigned to and N^, respectively. 

To obtain a frequency response, various values are assumed and assigned to 
XgCA , 00 ) , a vector that is a function of both frequency and amplitude . Then the 

vector signal is traced for each calculation . The phase and amplitude of each 
successive vector are noted through each block and summation junction until an 
input-output vector relationship is determined. 

The following logic was programed on a digital computer: 

s = jco , where 00 is variable ^6 ^ 

Xg variable 

^5 = ^V2^5^2 
X^-X^s<N^ 


7 s 




^9 = '^ 1 ^ 6^2 


P 
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APPENDIX E - Concluded 


X 


This makes g-^ a function of and co . The quantity Xg can be any one of the 

!p 

parameters 0, 0, *0, or a^, depending on the transfer functions Tg, Tg , and 
where 


N. 


Nr 


rate limit, 15 degrees per second or 12.6 degrees per second 
position limit , 2,5° or 6.5° 

705. 6(s^ + 67.8s + 2553.5) 


actuators , 


(s^ + 50.5s + 1568)(s + 33.9)= 


0 


■6.00(s + 0.8) 


^ s^ + 1.50s + 4.00 

0.375(s + 8) 


shaping , 


s + 4 


For a /6 at the center of gravity, 

YX 6 


T3-0 


T4«0 


T - Y V- 4.5(s^ + 0.60s - 27) 
® ^ s^ + 1.50s + 4.00 


For 0/6 at the cockpit, 

6 


r=^B= 


3 6_ ^2 


e s^ + 1.57s + 246 

T =i- 
4 s 


Te = o 
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APPENDIX F.-FEEL SYSTEM ANALYSIS 


The feel system of a fighter airplane is usually represented by a second-order 
system (ref. 10) . However, flight data were available at the PIO frequency that 
indicated that there was a hysteresis loop in the function . It seemed reasonable to 
assume that the pilot applied a force to provide the desired stick displacement; there- 
fore , stick position, 6^ , was considered to be the output. Figure 23 presents data 

s 


-20 -16 -12 -8 -4 0 4 8 12 16 20 



Forward N 

Figure 23. Stick position versus stick force . Flight data; f=0.5cps; 

A5 

e 

-^^= 0.0863 deg/N (0.384 deg/lb), 
s 

typical of stick position as a function of stick force. The hysteresis loop shown by 
the straight lines was considered to be representative of the feel system. The 
breakout force was approximately 22.24 newtons (5 pounds) . The gradient 
A6 

e 

was approximately 0.0863 degree per newton (0.384 degree per pound) . 
s 

Figure 24 shows the describing function for hysteresis adapted from reference 11-. 
With this describing function and faired data from figure 23, the feel system was 
represented as a function of stick output, 6 . To account for gearing change with 

s 

stick position, the following cubic equation was used: 

5g = f/Sg \ = 0.002786^ ^ + 0.45566 (FI) 

p V s/ ®s 
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APPENDIX F - Continued 



Figure 24. Hysteresis-describing function (adapted from ref. 11). 

Figure 25 compares the results of this equation with experimental data from a 
preflight control cycle . The data show some hysteresis due to the actuator . This 
hysteresis was accounted for in the digital computer program's description of the 



Figure 25. Mechanical gearing. Preflight control cycle (SAS off) 
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APPENDIX F - Concluded 


airframe plus augmentation system. The cascade combination of feel system plus 
the gearing is shown in figure 26. The amplitude and phase are presented in terms 
of output amplitude. This is actually the input command, 8^ , from the pilot to the 

P 

actuators (fig. 26) . As expected at low amplitudes , the phase difference is large 

5 

e 

and the amplitude ratio of is low . With increasing amplitudes of 6 , the phase 

s ®p 

difference decreases and the amplitude ratio increases . 



0 .1 .2 .3 .4 


6^ , rad 
e 

P 



Amplitude and phase characteristics of the feel system. 
= 0.0015 rad/N (0.0067 rad/lb) . 
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